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Abstract 


TeV or Super TeV neutrinos are expected to originate at the Gamma ray burster ( GRB ) events 
in the universe. These neutrinos are expected to be produced from the photo-meson interaction of 
the protons in the GRB environment. In the usual picture, the protons in the GRB hreball undergo 
Fermi acceleration in the ambient magnetic field of the GRB to very high energy. These protons then 
interact with the MeV, KeV or meV photons to produce Delta particles those subsequently decay to 
produce very high energy neutrinos. In this note we focus on the 100 TeV neutrinos produced in a 
GRB in the fireball model. We try to estimate the production point of these multi TeV neutrinos 
from the center. The strength of the Magnetic held there. Lastly but importantly modihcation in the 
spectrum because of neutrino absorption in the ambient magnetic held. The strength of the ambient 
magnetic held may change the neutrino spectrum. 


1 Introduction 

It is established now that GRB’s are among the most violent energy releasing events of the universe. 
Total energy release in these events could be equivalent to few solar mass and are seen to be distributed 
isotropically all around the universe. They are very high redshift objects occurring in nature - implying 
their cosmological origin. Since their discovery, the physics of the process that powers these events is yet 
to be understood. For some time now it has been conjectured that GRBs are source of very high energy 
cosmic rays. Moreover they have also been tipped off as one of the source of very high energy neutrino. 
The energy reachable in the GRB mediated events are several orders of magnitude higher than those 
available in the laboratory, that is one of the motivations why there is some interest in the community 
to study the physics of GRB. 

The origin of high energy particles in a GRB is ascribed to Fermi mechanism in the ambient magnetic 
field in the GRB enviornment. As we discuss later in this paper the same magnetic field might be respon¬ 
sible for modifying the v spectra. This requires good understanding of GRB physics as well as, careful 
evaluation of GRB parameters like temperature, Magnetic held, jet Lorentz Boost factor, matter density 
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in the rest frame of the jet etc. We discuss them below one after the other. 

As already mentioned, GRB’s are very high redshift objects, implying, they must have happened during 
the time of formation of the first generation stars of our universe. They have certain distinct features, 
such as, the rapid variability time (10“^ Sec), the non-thermal nature of the spectra, the hardness of 
the spectra (~ 100 MeV ). The existence of large energy in the power law tails of the GRB spectra— 
above threshold gives rise to whats usually referred in the literature as compactness problem [1]. The 
naive estimate of the source size from the variability of the spectra runs in to problem. This is princi¬ 
pally because of the fact, that a stationary compact object with an energy release, of the order of 10^^ 
erg/sec would produce an optically thick environment; however the fact that one observes non-thermal 
spectra, gave rise to what is called Compactness Problem Of GRB. The time variability and hardness of 
the spectra together suggest the source to be extremely compact in size ( To ~ 10^ cm ), from which a 
highly energetic jet moves outward around the rotational axis of the object.Actually the unusually hard 
and non-thermal nature of the photon spectra suggest that the jet must be moving towards the observer 
with high Lorentz boost T, whose magnitude can vary between several hundreds to few tens of thousands. 

Inside the moving jet, the medium has appropriate optical thickness to last the radiation in the jet up 
to the variability time ( of the jet). That could vary between 1 to 10“'^ sec. To have the jet optically 
thick to ( One to two MeV photons ), it has to have some baryonic contamination in - composed of 
p and n, [2]. They are accelerated to very high energy in the shock and ambient magnetic held of the 
GRB through Fermi mechanism and subsequently would produce pions through photo-meson or internal 
collision processes. These pions hnally decay to high energy neutrinos, electrons and positrons. In this 
note we would be focusing mainly on neutrinos. 

The neutrinos produced in GRB, are emitted at various stages of the jet evolution with varying en¬ 
ergy. Starting from low energy neutrinos at the beginning to extremely high energy neutrino energy at 
the end. Typically GeV neutrinos are expected to be produced at the initial stage, followed by 5 TeV or 
more and then > 100 TeV ones. 

Neutrino properties are known to be sensitive to ambient Electro-magnetic held and matter held [3]. In 
particular there are several decay channels, forbidden in vacuum but are allowed in a background mag¬ 
netic held. Actually these are the motivating reasons to look into aspects of the magnetic held structure, 
medium properties; at locations where neutrinos of specihc energy would originate. With this purpose 
in mind, in this note we would estimate the radial distance form the center of the star {tt:) where the 
energetic pions form to decay into 100 TeV neutrinos. We then go ahead to estimate the magnetic held 
there, followed by estimates of magnetic held induced neutrino decay. 


1.1 Neutrino Production In GRB. 

As has been mentioned in the beginning, in this note we investigate ultra high energy neutrino production 
and their interactions in the ambient medium and the possibility of change in the spectrum because of 
that. The reaction dynamics is as follows, pions hrst decay to muon and muon type neutrino, i.e., 

7r+ ^ -I-(1.1) 

followed by the muons decaying to, 

^ e+ -b (1-2) 

From Eq.(l.l) and (1.2) its clear that, emission takes place from two different reactions. The 

hrst emission would be taking place within a radius ( the pionosphere radius) within which the pions 
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produced from photo meson interaction would decay. The other region being the location where the 
energetic muons decay (usually referred in literature as r^). For a given rd one can derive the energy of 
the muon anti-neutrinos observed on earth ( ) or vice versa. The muon anti-neutrino in Eq. (1.2) 

is an end product of a three body decay, therefore the energy carried by the in the rest frame of 
the muon- would be (roughly!). Where is the rest mass of the muon. Hence the anti-neutrino 
energy as observed on earth would be given by, 

epM = In units of )l=c=l. (1.3) 

Where F;, is the Lorentz boost of the wind with respect to an observer on earth. On the other hand, if 
the muons with half life originate in the jet wind of the GRB ( extending up to distance rd), then it 
follows from there that, 


d — FpF/it^. 


(1.4) 


From relations Eq.(1.3) and (1.4) one can relate the observed anti-neutrino energy with that of the wind 
dynamical radius, i.e.. 


Ob ^ Tn^rd 

■ 


(1.5) 


For the muon type neutrinos, produced in a frame where the pion is at rest, the available kinetic energy 
is roughly close to (m,r ~ w^) = 34 MeV and the resulting muon is non-relativistic. So even if most of 
the energy is carried by the it would not be very different from the or the Ve, i.e ~ Although 

their point of origin in the GRB environment would vary ( because of difference in lifetime of pion or 
muon), hence change in the Lorentz Boost, however given the detector efficiency, that probably would 
not be detectable. Therefore we would consider the point of their origin to be same. 

We would like to note, in absence of magnetic field the ratio of number density of neutrinos (uo^riofi) 
would be 1:1, at the production point. However in a magnetic background new decay channels, e.g., 
v —> opens up with nonzero probability a. In general a depends of lepton mass mi their energy 

Eo, strength of Magnetic field B etc. Since these many of these quantities are flavor dependent, implying 
different absorption rate for different neutrino flavors. And to find out the absorption rates we need 
to estimate the parameters like Strength of Magnetic Field at a particular location. The energy of the 
neutrinos being produced there etc. In view of this this has been estimated in this paper, following the 
fireball model of [2]. 

To estimate various physical parameters of a GRB, we would stick to the fireball model of GRB ( 
successful) and try to refrain from making any explicit assumption about the progenitor. The postulates 
and assumptions we make, would be based mostly on [2], [7], [5] and [6]. 

The organization of this document is as follows, in section B, following fireball model we discuss the 
physics of GRB. Section IB is devoted to estimating various parameters, and lastly section IV would 
deal with details of absorption rate and the asymmetry in number density of and for some specific 
energy > 100 TeV. we conclude by discussion and some numerical estimates of asymmetry. 


2 Physics Of GRB 

Numerical simulations based studies reveal that, the jet produced at the central engine of the GRB moves 
along the rotation axis with a Lorentz factor (Lf) F. After launch.jet consists of shocked jet material 
and shocked stellar materials. Initially, the jet heats up the region where it interacts with the stellar 
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material, till the pressure in this zone becomes comparable to the same in the jet. Following notations of 
[4] region ahead of shock frame is called up stream, the one behind the shock called down stream. As the 
pressure in the up stream termination shock grows material flows into the jet creating a reverse shock ( 
Lorentz factor F/j). Subsequently material in the shocked region down stream (jet head) move out ward 
with Lorentz factor Th- The synchrotron and inverse Compton is believed to diffuse in, thermals and 
power the movement of the jet head. The isotropic luminosity of the GRB usually is related with the 
temperature of the medium in the jet head.In the rest frame of the jet ( of Lf, F;i ) the retarded jet moves 
with a Lorentz factor [7]. 

After injection, during its passage through C/O/He core jet head moves with subreltivistic velocity. 
However beyond He core there is a significant drop in the density of the stellar material (envelope)i.e., 
p < 10“^ gm/cc. and the jet can accelerate. As will be discussed below, during this phase F/i oc jr, 
where r and Tq are the radial distance at the point of concern and the injection point. The jet accelerates 
till it’s L.f reaches a saturation value. Further out it follows a self-similar motion and decelerates while 
interacting with the (ISM);we note here that baryonic components (mainly protons) of the jet are accel¬ 
erated to high energy in the reverse shock and these high energy protons interacting with the thermal 
photons produce neutrinos. 

Near the outer envelope, jet is believed to be dissipation-less hence the jet fluid in GRB can be collision¬ 
less. It behaves like an ideal Magneto Hydro Dynamic (MHD) fluid. As a result it follows that the 
entropy and the energy of the fluid component can be held constant. Gonservation of entropy and energy 
implies [9], 

r^r(r)roT.^(r) = constant. (2.1) 

and 

r'^T{r)roT{r)T^{r) = constant. (2.2) 

In the equations above, Vo stands for radius at the injection point and its usually taken as To = lO^cm. 
We would like to point out that, r(r) the position dependent Lorentz factor, from now on, would be 
replaced by F/j unless explicitly mentioned. From the estimates of Luminosity in terms of temperature 
one obtains, the estimate of the temperature at the injection point turns out to be, 

To = (2.3) 

With L52 = 10~^‘^Liso erg/sec and Tq = 10“^ro7 cm. Assuming, r(r)|r=ro = 0(1), it is possible to 
estimate the rhs of eqns. (2.1) and (2.2). They turn out to be, 

r'^T{r)roT^{r) ~ 2.1 x (2.4) 

and 

r^V{r)roT{r)T^{r) ~ 2.5 x lO^^LgsMeV. (2.5) 

As before unless explicitly mentioned we would suppress the explicit position dependence of temperature 
T(r). Relations (2.4), (2.5) imply, 


F/i oc r. (2.6) 

Assuming the presence of only magnetic field, ((eR), in the rest frame of the plasma.), using dissipation¬ 
less fluid equations, it is possible to derive an exact relation, 


eB 


= rijmp 
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Where, we assumed the medium to be composed of heaviest charged particle available in the plasma, 
namely protons . Their number per unit volume in jet rest frame is defined as nj. This is justified 
because the contribution from electrons would be far too less. Further more if we assume them to carry 
the isotropic energy ( luminosity ), L ergs/sec, then the number density of protons, at some distance r ( 
to an observer on earth ) from the center of the GRB ( moving with Lf, T^ ) towards earth would be [6], 

For a locally charge neutral plasma, Eq.(2.8) would also provide the number density of electrons. 


The physics of the shock front is described by Hugoniot equations.They balance various thermodynamic 
quantities like, energy, pressure, etc., across the discontinuity of the shock front. For instance equating 
the pressure between the forward and backward shock, (with density of material p in units of gm cm“^) 
it has been shown [7] that , 


F 


4 

h 


4p 


IbTrr^p 


(2.9) 


Following [8] near the outer envelope one can parametrize p = p* 
Using this parametrization, it turns out that, 


1 )' 


with p* < 10 7gm/cc. 


Th 


Lr 

IbTrp* {R — r)^ 


( 2 . 10 ) 


which matches with the suggested form F^ oc r, if one expands the denominator in powers or for 
r < R and retains the leading order term. 


3 Parameters 


3.1 Temperature Of the plasma 

In this section we would try to estimate the temperature, that the shocked plasma jet moves with. Since 
photons are believed [2] to be radiated from the shock heated plasma jet, moving with Lf F^j then all 
the radiated photons from the GRB would account for its total luminosity Liso. Using the emissivity- 
temperature relationship for pure photon gas it is possible to establish the temperature of the emitted 
radiation. The same under this assumption turns out to be. 


T 

J. 



(3.1) 


However this relation is bound to get multiplied by correction factors (< 1) if the photons constitute 
a part of the total luminosity L. However to keep the physics simple we assume the photons to carry the 
total energy of the GRB. We conclude this section by noting that, the number density of the photons at 
this temperature would be given by. 


Uj = T^ (3.2) 

which essentially follows the laws of pure photon gas. One can use eqn. (3.1) to express photon number 
density in terms of total luminosity L and radial depth from the center. 
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3.2 Ttt: Pion Formation Radius 


In previous sections we have estimated the various parameters related to the production of high energy 
protons as well as the approximate temperature profile of the shocked plasma head. In view of these 
estimates it is appropriate to estimate the the distance from the origin where these protons ( of energy 
Cp ) would scatter off the thermal photons to undergo A formation. In the jet head, the protons at the 
resonant energy would scatter off thermal photons to produce A as long as,at the rest frame of the jet, 
the following energy condition, in obvious notation, is given by, 

e^ep = 0.2GeV^ (3.3) 

Before proceeding further we note that, the reaction rate R for p + ^ A would be given by 

R = n^njap^^^A [1 — cos (0)] ~ n^njap^^^A, (3-4) 

Here we have assumed 9 = 90°. Its appropriate to note that, in order to have resonant photon proton 
scattering to A, inside a GRB shock cavity, number of protons (rij) in unit volume should be equal to 
R X 6t; with 6t being variability time for the GRB, this leads to the condition. 


Th^(Tp ^—— 1 . 


(3.5) 


In eqn.(3.5 ) above, ap^^A = 5 x 10 ^®cm^, is the resonant proton photon cross section for A production. 
Eqn. (3.2) implies at a distance r from the center, rij = where T-y = e-y, i.e., the energy of the 

h 

photon. Using this relation, eqn. (3.5) implies, 

Liso O'p'y^A 1 /r, 

4^ p2r2 g = 1- (3-6) 

In order to have resonance scattering, the energy condition, i.e., eqn. (3.3) should be satisfied. Replacing 
in terms of Cp by using, (eqn. (3.3) ) and the expression for from eqn. (2.9), one arrives at. 


r 


4 _ 
h 


y i^R^LisoP^^ CTp.y —fA Cp 

v7 X (.2 GeV^) 


2 

5 


(3.7) 


If we assume, the variability time 6t ~ 10“^ sec, L^so = 10^^L52 erg/sec, and the energy of the proton 
tp = 10® GeV. Then the radial distance from the center to have pion formation, r^r turns out to be, 
~ 10^° ® cm for R = lO^^cm and = 10“^ gm/cc. For more energetic protons, the radius of pion 
formation would shift further out-wards. For instance production of neutrinos of energy > lO^^eV, 
would form further out. Once we have the estimate of the pion formation radius, we can estimate the 
strength of the magnetic held and hnd out the held induced modihcation to the spectra. The estimation 
of the magnetic held is done in the next subsection. 


3.3 Magnetic Field Estimation 

In this note we have assumed the system to follow ideal MHD conditions. That is to say that the magnetic 
energy density can be expressed following eqn. (2.4) 

~ 2.1 • 10®®i5/^ (3.8) 

Where we have used the energy density at the injection point i.e., r = To = 10^ cm with magnetic 
equipartition taken to be of order unity. We note that Tq is also the width of the shell, that the GRB 

^We have used the fact Liso^t = 10®^L52erg = 0.6 X 10®®L52 GeV 
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wind expands with. Conservation of magnetic field energy at the wind rest frame implies at any radial 
distance r, in the GRB environment, 

= dTTr^B^. (3.9) 

leading to the relation, B x r = Bo x ro implying B = Bo [roir). Using these relations, the magnetic 
field strength, in terms of GRB luminosity turns out to be. Bo = 4.2.10^^L52^]^/^ . Hence the magnetic 
field strength at distance r from the center turns out to be, 

H ^ . (3.10) 

For conditions prevailing at the pion formation radius, corresponding to 10^ GeV protons, the strength of 
the field- compared with the critical Field strength Bc{sim'ml)- i.e ~ 0.01. The numbers mentioned 
here can get an order of magnitude variation depending on the details. Since the purpose of this note is 
to get an order of magnitude estimate of the effects, we reserve the details for future communications. 


3.4 Estimates Of Energy Loss. 


The high energy from the decay of A would undergo energy loss by inverse Gompton (IC) 

scattering, if the optical depth is large. However even at low optical depth the end products remain 


coupled to photons provided the photon drag time to 


is comparable to their decay time [13]. 


One can verify that tn is indeed comparable to the decay time, ensuring that the extremely relativistic 
pions and muons would undergo IC losses during their life time in the jet. However during the same 
period, they are expected to be accelerated in the reverse shocked plasma, to gain energy. Recalling 
the expression of pressure at the shock front to be Pr = 4 rijmp, the equation of motion for the 

charged particles at the shock front would be given by: 


±=4 

dt 


rijirij X Area. 


If the jet makes an angle 9, at a distance r from the center then the area would be 7r(r0/2)^. 
further use the fact that one gets. 


One can 


Trrijmjr^ 


Trnjrl 


While coming to the relation above, we have made a conservative estimate, by replacing r by Tq and 
neglecting the rest mass of the particles. During the life time t, the fraction of energy gained, by the 
particles would be, 

AP, = (3.13) 

The inverse Compton energy loss, in the ultra relativistic regime has already been estimated in [12], and 
is given by 


A/c = 


ifj r ctt (we / rrij ) ^ n.y 


Finally the ratio ,^ This ratio probably would decide whether the IC loss is significant 
enough to cause any significant suppression in neutrino production. For the kind of situation we are 
interested in, it seems this ratio would be less than unity hence, not much of suppression in the high 
energy neutrino spectrum. However this is a very tentative statement and one still needs to perform a 
quantitative estimation. 
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4 Discussion and Outlook. 


High energy neutrinos moving in a magnetic field can decay to v ^ . Where I can be either electron 

or muon or tau. The important point here is to note that, such processes are otherwise forbidden in 
vacuum. It has also been noted that, these processes are insensitive to neutrino mass and mixing in the 
leptonic sector. In a purely magnetic field the probablity of these processes depends on the parameter, 

jy _ eBp±_ ^ Depending on whether V is greater than unity or less than unity, the probablity for this 

. . . . 

reaction varies significantly. In a magnetized vacuum the probablity for such processes, a (for V » 1), 
is given by 


Oil 


Gp {9I + 9a) 


InV - i|n3 


1 9 33 71^ 1%^ - 63^^ 
75§56r4[|] gl + g\ 


(4.1) 


We would like to remind here that 7^; is Euler’s gamma function and r[|] is the Gamma function. The 
last term in eqn. (4.1) contributes significantly when the neutrino flavor index and lepton flavor indices 
are different, we however would be interested in such a case. 


In a GRB, the location of the neutrino production of a particular flavor strictly speaking would be 
different. However since the muon life time is pretty much small, therefore the location of v production 
from 7r“'" decay and the same from decay though would be different however the difference in length may 
not create significant changes in the medium properties. Therefore, for all practical purposes, we would 
assume them to originate from the same location. If the number density of u of certain energy Ei, at the 
production point is given by (E) then after travelling a distance dx the modified number density would 
be given by dn^^{E) = —no,y^{E)adx. And the solution to the same eqn., gives n^^(^E) = nQi,^{E)e~°‘'^] 
where ai is the flavor dependent absorption coefficient. We can further quantify the flavor dependent 
absorption as. 


Ki = —In 


T-Vi(E) 


Therefore one can define the asymmetry parameter Ap, to be: 


Ap — 


Eu.. - 




av 


(4.2) 


(4.3) 


A measure of this quantity Ap can in principle give us information about the nature of the medium at the 
production site of the neutrinos as well as the neutrino properties. In section [3.2] we had estimated the 
pion formation radius for Cp > lO^GeV followed by the estimation of magnetic field at that location. It 
should be noted that any enhancement of the magnetic fields due to charge particle streaming has been 
ignored here. We assumed the magnetic field to be toroidal and the neutrino propagation direction to 
be perpendicular to it. Furthermore we assumed the strength of the magnetic field to be held constant 
over the thickness of the shell (in its rest frame) i,e Tq = lO^cm. This can be justified because any 
variation to magnetic field over a distance larger than ^ would not make any significant contribution 
to the probablity amplitude we are interested in. 


Coming back to the question of asymmetry, as has already been discussed the asymmetry parameter Ap as 
has already been discussed, is far from zero. For very high energy neutrinos it seems to decrease compared 
to the low energy neutrinos, see fig. [4]. However estimation of actual degradation of energy spectrum 
of TeV or Super TeV Neutrinos would require more careful estimation of the parameters discussed here 
as well as taking matter effects into consideration. However there is an interesting consequence of this 
effect, that is, ther might be some enhancement in the number of neutrinos coming from the decay of 
muons originating from the reaction Vf, —>■ gT. That might have interesting consequences for the 
upcoming neutrino experiments. 




Figure 1: Asymmetry parameter {Ap, see text) for Vg and Vp, for the High Energy neutrinos in the GRB 
enviornment. x - Axis corresponds to neutrino energy in MeV and y Axis corresponds to Ap. Magnetic 
field strength in units of electron mass is taken to be 0.001. Energy is taken to be same for both the 
flavors. 
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